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Recent experiments suggest that polarized photons may couple significantly to pseudoscalar par- 
ticles such as axions. We study the possible observational signatures of axion-photon coupling for 
radiation from magnetic stars, with particular focus on neutron stars. We present general methods 
for calculating the axion-photon conversion probability during propagation through a varying mag- 
netized vacuum as well as across an inhomogeneous atmosphere. Partial axion-photon conversion 
may take place in the vacuum region outside the neutron star. Strong axion-photon mixing occurs 
due to a resonance in the atmosphere, and depending on the axion coupling strength and other 
parameters, significant axion-photon conversion can take place at the resonance. Such conversions 
may produce observable effects on the radiation spectra and polarization signals from the star. We 
also apply our results to axion-photon propagation in the Sun and in magnetic white dwarfs. We 
find that there is no appreciable conversion of solar axions to photons during the propagation. 



I. INTRODUCTION 

The axion is a hypothesized pseudoscalar particle aris- 
ing from the breaking of a U{1) Peccei-Quinn symme- 
try, introduced to explain the absence of strong CP 
violationP, H i]. The axion is also an ideal candidate 
for cold dark matter that makes up five-sixths of all mat- 
ter in the universe. The allowed axion mass is in the 
range of lO'^ <ma< lO^^ gv[3|. 

A general property of the axion is that it can couple 
to two photons (real or virtual) via the interaction 

La-y^.^-^gaF^^pf^'' ^ga^-B, (1) 

where a is the axion field, F^i, [F^^") is the (dual) elec- 
tromagnetic field strength tensor, and g is the photon- 
axion coupling constant. Accordingly, in the presence 
of a magnetic field B, a photon may oscillate into an 
axion and vice versa. Exploiting such photon-axion os- 
cillation, experiments based on "photon regeneration" 
("invisible light shining through walls") set the con- 
straint g < 6.7 X 10"'^ GeV"^ (95% CL) for axions with 
■nia < 10^3 evd). The first results from the CERN Ax- 
ion Solar Telescope (CAST) experiment imply an upper 
limit g < 10~^° GeV"^ for < 0.02 eV 'g. Other ex- 
perimental and astrophysical constraints on the axion pa- 
rameters (mass 771(2 and g) are reviewed in Refs. 0, 0,ll|- 
Recently, the PVLAS collaboration has reported mea- 
surements of the rotation of the polarization of photons in 
passing through a vacuum cavity in a magnetic field [9|. 
If interpreted as due to the conversion of photons to ax- 
ions, this may imply g ~ (1.6 — 5) x 10~^ GeV~^ for 
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nia ~ lO"'^ eV. Unfortunately, this result contradicts the 
constraint from the CAST experiment and the previous 
energy-loss limit from globular cluster stars. Thus, ei- 
ther there remain systematic effects in the PVLAS ex- 
periment or nonminimal models of psuedoscalar-photon 
coupling are required IQ}- 

In this paper, we study the effects of axion-photon 
coupling on radiation from magnetic stars (see [111 ] for 
previous works). We will consider generic psuedoscalar 
parameters with nia ^ 10"'^ eV and g ^ 10~^ GeV~^, 
not necessarily the QCD axion (for which there exists a 
unique relation between nia and g). Our goal is to un- 
derstand and quantify what might be the observational 
signatures of axion-photon coupling for various axion pa- 
rameters. Our method and basic results (such as the 
axion-photon conversion probability) can be applied for 
all axion parameters and in different astrophysical en- 
vironments. Most of the paper will deal with magnetic 
neutron stars, for which most significant effects are pos- 
sible, but we also discuss applications of our results in 
the contexts of magnetic white dwarfs and the Sun. 

Our paper is organized as follows. In Sect. II we study 
axion-photon propagation in the vacuum region outside 
a magnetic neutron star, taking into account of the varia- 
tion of the magnetic field. Section III examines the prop- 
agation in a magnetized plasma characteristic of a neu- 
tron star atmosphere. We study in detail the properties 
of the axion-photon resonance as the photon-axion prop- 
agates in an inhomogeneous medium. Depending on the 
axion coupling strength and other parameters, significant 
axion-photon conversion may take place and such con- 
version may leave an imprint on the radiation spectrum 
and polarization from the star. In Sect. IV we consider 
axion-photon propagation in the Sun and in magnetic 
white dwarfs, and Sect. V contains a brief discussion of 
our results. 
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II. AXION-PHOTON PROPAGATION IN AN 
INHOMOGENEOUS MAGNETIZED VACUUM 

Here we consider the effect of axion-plioton coupling 
as tlie plioton propagates radially outwards through the 
vacuum region outside a magnetized neutron star. The 
stellar magnetic field is assumed to be dipolar. Both of 
these assumptions are valid several radii away from the 
star. 



where = ma/(10~^ eV) and gg = g/{10^^ GeV^^). 

Note that Eq. ([2]) is valid in the weak-dispersion limit, 
i.e., A||/w, |Aa|/w and Am/^^ are all much less than 
unity. 

In the absence of photon-axion coupling (Am = 0), 
and with -E||,a oc e"**^^, we find k/oj = 1 + A^^/ui and 
1 + Aa/uj for photon (the parallel-mode) and axion, re- 
spectively. 



A. Basic Equations 

In an external static magnetic field B, the evolution 
of the photon electric field E and the axion field a of a 
given energy iv (so that E, a oc e*'^*) takes the form [ll| 
(in units with ?i = c = 1): 



dz 



a 

Eh 



UJ + Aa 

Am 



^ Ail 



a 

Eh 



(2) 



where E\\ is the photon electric field in the plane spanned 
by B and the z-axis, and we have assumed that the wave 
propagates in the z-direction. The matrix element A|| 
arises from vacuum polarization [13, [13, 13, [3] and is 
given by (for lo <C me, the electron rest mass) 



A|| = — quj sin^ 



(3) 



where 9 is the angle between the direction of propagation 
and the magnetic field, and g is a dimensionless function 
of 6 = B/Bq, with Bq = mlc^/{en) = 4.414x lO^^ G the 
critical QED field strength. For 6 < 1, q = 7ab'^ / {4:5tt). 
A general fitting formula (accurate to within 3% for all 
&'s) which has the correct 6 <C 1 and 6^1 limits is [l^ 

1.2b 



7a 2 - . , - 
q = -r^o q, with q 



1 



1 + 1.336 + 0.5662- 



(4) 



Thus 
A 



■II 



7a .^2 



0.1807q/3^wi eV = 9173g/3^wi cm-\ 

(5) 

where uoi ^ uj/{l keV) and 



B 

d ^ b sin 6 = sin I 

Bo 



(6) 



The matrix element Aa in Eq. ([2]) is due to the finite 
axion mass, and Am arises from the photon-axion cou- 
pling: 



A„ 



= -5 X lO^^'^m^cjf^eV 

= -2.538 X lO^^m^wf ^cm"\ (7) 

AM^^gBsine^^f] 

= 1.306 X 10"'^(79/3 eV 

^ 6.627 X 10^^59/3 cm-\ (8) 



B. Evolution in a Varying Magnetic Field 

It is convenient to introduce the quantity Afc and the 
mixing angle Om via 



2 ll/2 



and 



Afc = 2[(Aa- A||)V4 + A2,] 



tan 20m = T-. 1 — :-r- = A. 



(Aa-A||)/2 

Then the evolution equation Q can be written as 



(9) 



(10) 



dz 



a 

Afc / cos26i„ 
2 \ sin26'„ 



All + Aa 



sin 2dm 
- cos 29 n 



ill) 



where I is the unit 2x2 matrix. 

The eigenvalues and eigenvectors of Eq. pT|) are 



Aa Afc 



2 2 
for the mode, and 



a 

Eh 



Aa 



Afc 



a 

En 



cost 
sin^ 



— smfy 
cos 9„ 



(12) 



(13) 

for the "— " mode. 

Figure [T] shows an example of the mode eigenvalue as 
a function of (3 (with the other parameters fixed). For 
sufficiently large /3 (such that A|| ^ Am, |Aa|), we have 
Afc ~ A|| — Aa, thus fc+ = + A||, fc_ ~ w Aa and 
tan 29m — 0. As f3 decreases, | tan 29m \ first increases and 
then decreases. For sufficiently small /? (such that | Aa| ^ 
Aj\/, A||), we again have Afc ~ A|| — Aa, fc+ = w + A||, 
fc_ ~ w -I- Aa and tan26'm ~ 0. Thus, for both large 
and small /?, the mode always represents the axion 
and the "— " mode always represents the photon. Strong 
mixing may occur for intermediate /3's when tan 29m ^ 1, 
but there is no "level crossing" . 

Consider a photon created in the large-/? region of a 
neutron star. As the photon propagates outward to the 
low-/3 region, it may partially convert to the axion. To 
calculate the net conversion probability, we define the 
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FIG. 1: The solid curve shows Afc/2, and the dashed curve 
shows \6'^\ (assuming that the field is dipolar and that at 
the stellar surface (3 = Ps — 1 and R = 10^ cm). The light 
dot-dashed line shows Am and the long-dashed line shows 
(A[| — Aa)/2. The parameters are: nia = 10~^ eV, g = 
2 X 10~^ GeV~\ Lo — 1 keV. The "Strong mixing" region 
(between the two vertical solid lines) is where (Ay — Aa)/2 < 
Am or P- < P < P+; nonadiabatic mode evolution occurs 
in the region where Afc/2 < \9'^\ (cross-hatched area). The 
vertical dashed line denotes /3* where the maximum mixing 
angle is achieved. 



mode amplitudes to be and and expand the 
photon-axion system in terms of the normal modes: 



= A, 



cos 6„ 
sinf?„ 



— sin On 
cos Om 



(14) 



Submitting the above in Eq. pip , we find that the mode 
amplitudes satisfy the equation 



d 



'dz \ A 



A+\^( AA:/2 \( A 



(15) 



where ' denotes derivative with respect to z. In deriv- 
ing Eq. (|15p we have neglected the non-essential diagonal 
term (proportional to I) in Eq. pT|) . 

In deriving the evolution equation (Eq. llSp we have 
assumed that the orientation of plane spanning the pho- 
ton wave number k and the magnetic field B does not 
change along the ray path. This is valid for a static 
dipole field far away from the star (more than a few stel- 
lar radii away) and well inside the light cylinder radius 
(c/fig, where Vis is the rotation rate of the star). For a 
rotating star, the orientation of the k-B plane changes 
along the ray path, but this is a rather small effect if the 
rotation rate is much less than the maximum (breakup) 
value. Note that the static dipole approximation breaks 
down for r beyond the light cyclindcr radius. 

The mode evolution will be adiabatic if the condition 



7 ' 



Afc/2 



> 1 



(16) 



is satisfied. From Eq. PTI)1 . we find 

sin4^A' A' P' Ag + A,, 

where in the A' /A expression we have set g' ~ for 
simplicity. For a dipolar field, we have 



(18) 



where /3s is the field strength at the stellar surface, 
R = 10^i?6 cm is the neutron star radius. Thus (3' /P ~ 
-3/r = -(3/i?)(/3//3s)^/^. If the adiabatic condition 1^ 
is satisfied at every point along the photon's path, then 
the photon will traverse the region of varying B without 
any net conversion to the axion; otherwise, the photon 
may experience conversion and the emergent photon in- 
tensity will be smaller than the emitted intensity. 

It is important to map out the parameter regimes for 
which net nonadiabatic photon-axion conversion may oc- 
cur. To this end, we note to have any appreciable con- 
version, the photon must pass through a "strong mixing" 
region, where | tan20m| is not much less than unity. In- 
deed we see from Eq. pT|) that if | tan 20^ | ^ 1 along 
the photon path, then E\\ oc exp(— i / dz fc+) and no con- 
version occurs regardless of whether the condition (jl6p is 
satisfied or not. We define the "strong mixing" region by 
|tan26'„| > 1, or /3_ < /? < /3+ (sec Fig. [J), where 



/3± = 



lauj 



l± Wl- 



(19) 



In the above equation, we have used q = 1 since for the 
parameters regime of interest, /3± ^ 1 is well satisfied. 
Thus the "strong mixing" region exists only when 



9jK 
ma 

99 



> 



> 



rria 
99 

7715 



7a 
^457r 

0.728 



1/2 



(20) 



At the threshold point, /3± collapse to 



f3* 



777 



^ = 52.6^ 

Lu \ 7a LO 



5.26x10-^^, (21) 



and (Afc/2)* = m2/(V2w). 

In general (not only at the threshold point) the max- 
imum mixing angle occurs where /3 = /9* and has the 
value 



I tan 20™, Imax — — I 



99_ 
7715 



99_ 

m5 



(22) 



Suppose g/TUa > (g/rria)*, then the region where the 
adiabatic condition (|16p is most likely to be violated is 
around (3 — (3±, where \9m\ = n/8 and |sin46'm| = 1 
(see Fig. 1). First consider the the parameter regime 
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FIG. 2: The parameter domain for which nonadiabatic 
photon-axion conversion is possible. Above the sohd line 
[Eq. (|20|l ] (the unshaded regime), there exists a "strong mix- 
ing" zone (|tan2Sm| > 1) between photons and axions. In 
the unshaded region, we have gi/mi > (31 /mi),, and well 
inside the unshaded region the asymptotic results of Eqs. (|23p 
and (|24l) apply. In this (unshaded) region and below the 
long-dashed line, mode evolution is nonadiabatic at /3 = /3_ 
[Eq. (|24p ]: below the short-dashed line, the evolution also be- 
comes nonadiabatic at /3 = /3+ [Eq. (|23|) ] . The heavy lines are 
for / = Lolp-^R^^ = 1, and lighter lines for / = 10. 

g/rua ^ {g/iria)*- In this case, /?+ is determined by 
Am ~ A||/2 > |Aa|, which gives (3+ ~ 90TTgml/{7auj). 
The nonadiabaticity condition Afc/2 < |0^„| at /3 = /3+ 
then translates to 

59 < 0.755(^?/37ii?6-3)i/5^ (23) 

At /3 = /3_, which is determined by Am — |Aa/2| > 
|A|||, giving /3_ ~ ml/{2ujgml) = 1.9 x 10-^mlg:[^Lj^\ 
the nonadiabaticity condition Afc/2 < \9l^\ translates to 

g,<Um-Hu^lp-^R-'). (24) 

Figure [2] illustrates the parameter domain for which 
nonadiabatic photon-axion conversion may be possible. 



where we have assumed that /3 decreases monotonically 
along the photon path; this is appropriate more than a 
few stellar radii from the stellar surface. At large (3, we 
consider a pure photon state, with = 1 and = 0. 
We integrate Eq. ([25|) toward small /3 and obtain the 
asymptotic (/3 — > 0) , and the net conversion probabil- 
ity is then Fconv = [A+j^. 

Figures [3]|n] give several examples of our numerical in- 
tegrations. Figures consider the regime {g /ma) > 
(g/ma)*- At low energy (Fig. [5]), the mode does not 
evolve through a nonadiabatic zone [neither (|23p nor ((24|) 
are satisfied], thus the mode conversion probability is 
small. At intermediate energy (Fig. [3]), the condition 
([M)) is satisfied while is not, we obtain a large Pconv 
At large energy (Fig. |4]), both ([24]) and (|23|) are satisfied, 
and yet the net conversion probability is small. Fig. [J] 
gives some hints as to what is happening. When the wave 
travels through a non-adiabatic zone, the amplitudes of 
the modes change. There can be two nonadiabatic zones, 
one near /9_|_ and one near However, the sign of 9'^ 
is different in these two zones; therefore, the sign of the 
mixing term D in Eq. ([25)) differs as well. When both 
zones are active (as for w = 20 keV in Fig. 2]), some of 
the changes in the inner zone (/3+) are undone in the 
outer zone (/3_). With this in mind, the peak conversion 
as a function of energy lo should occur approximately 
where Eq. (^5]) holds but Eq. breaks down, so that 
the non-adiabatic region around /3_ is as large as possi- 
ble without the region around /3+ being non-adiabatic as 
well. Thus, for a given axion parameters g, rua satisfying 
Eq. (j20p , the maximum conversion occurs for the photon 
energy between 

CO ~ 0.3g9/'m2(/3,i?3)i/2 (27) 

and 

uj^2gl/\psRiy^^ keY. (28) 

This nonmonotonic behavior of Pconv is shown in Fig. [T] 
As noted before, for (g/ma) < {g/nia)*, no "strong mix- 
ing" region exists in the photon's path, the net conversion 
probability will be small, regardless of whether condition 
(fTB)) is satisfied or not (see Fig. [5]). 



Numerical Calculations of Photon- Axion 
Conversion Probability 



D. Effect on Radiation Spectrum and Polarization 
fron Magnetic Neutron Stars 



To obtain the photon-axion conversion probability, it 
is necessary to perform numerical calculations. To this 
end, we rewrite Eq. pS]) as 



dln/3 



A+ 
A. 



C iD 

-iD -C 



where 



C 



Afc dr 
~2~dln^' 



D 



A, 
A. 



rfln/3' 



(25) 



(26) 



To illustrate the effect of photon-axion coupling on the 
radiation from magnetized neutron stars (NSs), we con- 
sider an example of a NS with surface magnetic field 
B = 10^^ G and surface teinperature = 0.5 keV (see 
Fig. [11). As shown in Refs. O, [3, due to the vacuum 
resonance effect [l^, [2l| occurring near the NS surface 
(but outside the photosphere) (see also Sect. IIIll below), 
the surface radiation from such a NS is dominated by the 
extraordinary mode (or _L mode) for ^ 1 — 2 keV and 
by the ordinary mode (|| mode) for w > 1 — 2 keV. In our 
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FIG. 3: The upper panel shows Afc/2 (solid line) and \0'^\ 
(dashed line), similar to Fig. 1. The lower panel shows the 
evolution of mode amplitude, starting from A- = 1, A+ = 
at a small radius (large /3). The parameters are rUa = 
10"'^ eV, p = 2 X 10"" GeV-\ oj = 1 keV, Ps = 1, and 
i?6 = 1. 



FIG. 5: Same as Fig. 3 except lo — 0.1 keV. 
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FIG. 4: Same as Fig. 3 except = 20 keV. 
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FIG. 6: Same as Fig. 3 except ma = 10 eV, g = 5 x 
10"" GeV-\ u = 1 keV, (is = 1, and = 1. In the 
upper panel, the light dot-dashed line shows Am and the 
long-dashed line shows (Ay — Aa)/2. Unlike the cases shown 
in Figs. ([3)l-([5|, no "strong mixing" region exists in this case. 
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FIG. 7: Photon-axion conversion probability Pconv in vac- 
uum as a function of photon energy. The parameters are 
rria = 10"^ eV, /3s — 1, Re — 1; the solid line is for g — 
2 X 10"^ GeV"\ the short-dashed line for g = 10"^ GeV"\ 
and the long-dashed line for g — 4 x 10^^ GeV~^. 

calculation leading to Fig. [U we assume (for simplicity) 
that the total intensity near the NS surface is g iven by 
a blackbody with temperature (see, e.g., for 
more accurate atmosphere model spectra), and we cal- 
culate the intensities of the two photon modes {I± and 
/|| ) outside the vacuum resonance layer using the method 
described in ^20] . With no photon-axion coupling, the ob- 
served photon intensity is I± + /|| , and Stokes parameter 
Q is I±— /|| (Note that the Stokes parameter depends on 
the orientation of the axis used in measuring linear polar- 
ization, but one can always choose the appropriate axis 
so that U = 0; Refs. give a more detailed descrip- 

tion of the phase-dependent polarization signals). With 
nonzero photon-axion coupling, as the photon propagates 
from the stellar surface to the observer, the ||-mode inten- 
sity will be reduced to /||^ = (1 - -Pconv)^||, where 
is the photon-axion conversion probability calculated be- 
fore (Fig. [7|), while the _L-mode intensity is unchanged. 
Thus the observed total radiation intensity is I± + I"^ , 
and the Stokes Q parameter is I± — . 

For NSs with Bs>7x 10^^ G, the radiation from the 
NS surface is dominated by the _L-mode, which is unaf- 
fected by any further photon-axion coupling as the pho- 
ton propagates to the observer. Thus, for such NSs, the 
axion would have a much smaller effect on the observed 
NS radiation than depicted in Fig. [H 



III. AXION-PHOTON PROPAGATION IN AN 
INHOMOGENEOUS PLASMA 

In this section we consider the effect of photon-axion 
coupling on the photon propagation in the atmospheric 




G) (keV) 

FIG. 8: Spectrum (upper panel) and Stokes parameter Q 
(lower panel) of a magnetized neutron star with surface mag- 
netic field Bs = 10^'^ G and temperature Ts = 0.5 keV. 
In both panels, the heavy lines show the results including 
photon-axion coupling with the parameters rua = 10~^ eV 
and g — 2x 10~^ GeV~^, while the light lines neglect photon- 
axion coupling. In the upper panel, the solid lines show the 
total radiation intensity, the short-dashed line shows the _L- 
mode intensity and the long-dashed line shows the ||-mode 
intensity. 

plasma of a magnetized NS. Here the magnetic field is 
constant, but the plasma density varies on the scale of 
centimeters (the atmosphere scale height at Tg ~ 10^ K). 

A. Vacuum resonance 

We first consider wave propagation in a magnetized 
plasma, including the vacuum polarization effect, but ne- 
glecting axion-photon couphng 18, 20, 2ll [H H [H]. 
In the weak-dispersion approximation (i.e., the index of 
refraction of the photon is very close to unity), the elec- 
tromagnetic wave equation in the plasma including vac- 
uum polarization (but no axion) is given by Eq.(30) of 
[r§ |. Using the sign convention adopted for this paper, 
we have 

,±( ^\\\-^( 2 + ^^11 <^i2 \( E\\\ (29) 

Neglecting the protons (i.e. setting the proton mass 
nip — > oo) and damping terms in the dielectric tensors 
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(more general expressions are given in [T8|). the matrix 
elements are given by 



(Til = (q- Ve) sin^ 9 - 
(722 = — TOsin^ 9 



1 



■ cos 



1 



Ue 



1712 



1/2 

-CT21 = «^ cos ( 



1 



Up. 



(30) 
(31) 

(32) 



Here Ue 



w^g/w^, with LOce = eB/{mec) — nieC^b the 



electron cyclotron (angular) frequency, and = Lo^f./u 

with Wpe = (47rneeVTOe)^/^ = 28.71 (YePi)^/^ cV the 
electron plasma (angular) frequency, where ne — Y^p/mp 
is the electron density, p is the mass density, is the 
electron fraction and pi = p/(l g cm~^). The quantities 
g, m are the vacuum polarization parameters: q is given 
by Eq. (HI) and [H 



m 



4a 
'45^ 



1 + 0.7265/4 + (4/i5)fo2 



(33) 



The weak-dispersion approximation adopted in Eq. (I29[) 
is valid valid for <C 1 and 6 ^ Sir/a ^ 10'^. 

With cx: e~*'^^, Eq. (|29p determines two normal 

modes. For ^ I, the modes are almost linearly polar- 
ized, except at the "vacuum resonance" [l^, [2l|, [23, [13, 
[25j |. which occurs at cth = CT22, or 



(9 



1 



(34) 



which for ^ 1, becomes = q + m. For a given 
photon energy, the resonance density is 

Pres(vacuum) = Q.964Y-\Bi40JiffB^ g cm'^, (35) 

where = B/{10^'^ G), uji ^ uj/{l keV), and fg = 
[{ab^ / 15tt) / (q+m)]^ is a slow varying function of and 
is of order unity (/s ~ 1 for 6 < 1 and /s < 4 for 6 < 100; 
we will use /s = 1 in remainder of this paper). At the 
vacuum resonance, the plasma and vacuum polarization 
effects are comparable, and both modes are circularly 
polarized. 

The evolution of photon modes around the vacuum 
resonance in an inhomogeneous NS atmosphere (varying 
density but constant magnetic field) has been studied 
before [H, . The key result is that away from the res- 
onance, the mode evolution is highly adiabatic. At the 
resonance, the mode evolution depends on the adiabatic- 
ity parameter 



ad) 



with 



Wad 



2.55 (/b tan6')2/3F-i/3 



(36) 



(37) 



where H — p/\p'\ is the density scale height along the ray 
(evaluated at the resonance point), and Hi = H/{1 cm). 



For an ionized hydrogen atmosphere, H = p/\p'\ — 
2kT/{mpg^ cos 6 kg) = 1.65 r6/(5*, 14 cos 6*^,3) cm, where 
T = lO^Tg is the temperature, = 10^'*(;*,i4 cm s~^ is 
the gravitational acceleration, and 9kg is the angle be- 
tween the ray and the NS surface normal. Since nona- 
diabaticity can only occur at the vacuum resonance, we 
can use the Landau-Zener formula to calculate the nona- 
diabatic jump probability between normal modes [iTl] : 



-L nimn — 



(38) 



Thus, for 7i.cs S> 1, the evolution is adiabatic, and the 
photon can convert from the ||-mode (ordinary mode) to 
the _L-mode (extraordinary mode) and vice versa. Var- 
ious implications of the vacuum resonance phenomenon 
for NS surface radiatioii spectrum and polarization are 
studied in Refs. [H, [M US IHI, HI • 



B. Axion-Photon Resonance 

We now consider the effect of photon-axion coupling 
on photon propagation in an inhomogeneous magnetized 
plasma. Combining Eqs. ^ and ([^5]) . we find that in the 
weak-dispersion approximation, the evolution equation 
for the photon (with two polarization components, 
and i?j^) and axion takes the form 




Am 

A|| 4 





CT12W/2 



(T2ia;/2 a;-t-CT22w/2 



where 




(39) 



(40) 



In Eq. Aq, A|| and Am are the same as in Eq. ©, 
CT22, cri2, (T21 are given by Eqs. (I3l|)-(l32]), and 



A„ 



UJV„ 



with F 



cos^9 
1 — u, 



+ sm^ 



(41) 



Obviously, in vaccum (ug = 0), Eq. ([M)) reduces to 
Eq. ([2|); with no photon-axion coupling (Am = 0), it 
reduces to Eq. (|29| . 

The axion-photon resonance occurs when the first and 
second diagonal matrix elements in Eq. p9p are equal, 
i.e., Aa = A|| + Ap, or 

qsin^ 9 + {ml/u'^) , . , , x 

Ve = (axion — photon resonance). 

^ (42) 
Equations are general, and are valid for differ- 

ent magnetic field regimes, such as those found in NSs, 
white dwarfs and the Sun. 

We now consider magnetic field strengths relevant to 
NSs and assume ^ 1. We consider a typical angle 9 
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that is not too close to 0° or 180°. For the parameters 
(axion mass, photon energy and field strength) of inter- 
est, gsin^6' '3> rnj^/u'^ (or A|| ^ |Aa|) is well satisfied. 
Then the axion-photon resonance is at 



(axion — photon resonance) 



(43) 



Note that the vacuum resonance is at We = 9 + m. The 
two resonances are well separated, and we can treat them 
independently — we justify this in Sect. IIII CI 

The vacuum resonance has been discussed in 
Sect. IIII XI We now consider the axion-photon resonance. 
Neglecting the component of the photon, the basic 
equation takes the form 



dz 



a 

En 



LJ + Aa 



A 



M 



LO - 



a 

En 



(44) 

Similar to Eqs. ([9]) and (fTO|) . we define Afc and the mixing 
angle via 



Afc = 2 [{Aa 
tan2^?„, — 



A 



A ^2 

M 



/4- 



1/2 



(Aa- 



A,)/2- 



The eigenvalues and eigenvectors of Eq. (|44p are 



fc± 



An+Aa + A^ 



± 



Afc 



cos ( 
sin t 



\I - 



- sm( 
cos 



(45) 
(46) 

(47) 
.(48) 



If we plot fc± as a function of the plasma density, we find 
the two eigenstates tend to "cross" at Aq = A|| -|- Ap, 
where Qm — ±45° — this is the axion-photon resonance. 
Similar to Sect. ITTBI we find that the mode evolution 
depends on the adiabaticity ratio 7 = Afc/(2|6''„J). As- 
suming that the variation comes only from density (while 
the magnetic field is constant), we have 



Ap p' 



" 4AAf p 
Thus the adiabaticity ratio is 
Afc/2 



sin^ 2B„ 



(49) 



7 



|Ap sin^26'„ 



4Al,i7 



1^ 2e„ 



(50) 



where in the last equality we have used A|| ^ |Aa|, and 
Pics IS the axion-photon resonance density (for a given 
photon energy lo): 



Prcs = 2.251; ^{BuUJifq g cm ■ 

At the axion-photon resonance p — prcs, \0n 
find 



AAl.H 



1.915 X 10" 



ujiq 



(51) 
7r/4, we 

(52) 



Note that 
sin 26„ 

with 



M 



1/2 ' 



(53) 



^3.614x10-^^1^. (54) 
Loq OLOiq 

For |1 - Velq\ > 4M/ sin^ 9, we have 



7~6.34x lO'^ 



59 



Ve/q 



(55) 



Thus, away from the resonance, the adiabatic condition 
(7 3> 1) is well satisfied, and nonadiabaticity can only oc- 
cur very close to the resonance. We can use the Landau- 
Zener formula to calculate the nonadiabatic jump prob- 
ability [13 



ump 



(56) 



Note that because of the level crossing, adiabatic mode 
evolution (71.cs ^ 1) corresponds to a conversion of the 
1 1 -mode photon to the axion and vice versa. Thus the 
axion-photon conversion probability across the resonance 

is ^conv — 1 -fjump- 

Comparing Eqs. (I36l)-(l37l) with Eq. it is inter- 

esting to note that high-energy photons tend to be more 
adiabatic at the vacuum resonance and but less adiabatic 
at the axion-photon resonance than low-energy photons. 

In Sect. IIII D I we shall present numerical calculations 
that confirm the analytical consideration given above. 



C. More Rigorous Treatment of Three-State 
Mixing 

Here we consider the eigenmodes of the three-state sys- 
tem, Eq. ([59)1 . We will show that the vacuum resonance 
and the axion-photon resonance can be considered inde- 
pendently of each other. 

The 3x3 matrix in Eq. ([55)1 can be written as (for 
We > 1) 



ujI + ujqG, 

where I is the unit 3x3 unit matrix, and 

AM 
M \{l-V)iiin^e -iV cose/{2uJ'^) 
_ iV cose/{2uJ'^) {2fh/lq)sm^e 



(57) 



(58) 



with 



V = Ve/q, 
A = Aa/{ujq) < 1, 

M = AM/iujq) 



(59) 
(60) 

3.614 X 10-^^1^ «1. (61) 
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FIG. 9: Eigenvalues of the axion-photon system in the 
plasma+ vacuum medium as a function of the plasma density 
parameter V = Ve/q, where Ve = [upejiSf' oc p. The insert 
shows the blowup of the region around the axion-photon reso- 



nance. The parameters are chosen so that A, M, u, 



^1/2 



< 1. 



Let (a, , £^j_) cx e 
via 



-ikz 



, and we define dimensionless K 



ujqK. 



(62) 



Obviously, K depends only on the dimensionless param- 
eters, y, A, M, Ue- Numerical values for K asa. function 
of V are shown in Fig. [9l 

— 1/2 

Since M <C 1, Ue <C 1, it is clear that away from 
the two resonances, the eigenvalues are simply 

K = A (axion), (63) 
K = {l/2){l-V)sin^ e (II -mode), (64) 
K = {2rn/lq) sin^ 6 (_L -mode). (65) 

We now consider the mode properties in the vicinity 
of the axion-photon resonance, |y — 1| <C 1. The third 
eigenvalue is simply K = (2m/7g) sin^ 0, describing the 
_L-mode of photon. The other two modes that "intersect" 
have eigenvalues that satisfy |i^±| <C 1, and are given by 



K± 



A-C l-V 



■ sm 



where 



C 



4 



(AMf 



I 1/2 



7? 



8rhup tan 9 



«1, 



(66) 



(67) 



and Vav is the value of V at which the axion-photon res- 

2{A + C) 



^ ap 

onance occurs 



sm 



(68) 



Note that if we set C = 0, Eq. §Eil reduces to Eq. (li^ . 
We see that effect of including the coupling with third 
mode (the _L-mode) is to shift the axion-photon reso- 
nance location by a small amount. The eigenvectors near 
the resonance are 



oc 



/ cos 






f - sin 9rn \ 


sin 




, CX 


COS Om 


V El 


+ / 







where 



E 



±± 



7qV COS0 

-i 777^ £-11 



\. 1/2 . 2a ""11=^' 
AmUe sm 

with \E^± \ < \E\\±\, and 



tan 29,ri — 



AM 



{V - Vap) sin 



(69) 



(70) 



(71) 



Again, the above expression reduces to Eq. (j46|) for C — 
0. We see that except for the slight shift of the resonance 
location, AA' = — K- and 0'^^ are the same as those 
given in Sect. IIII^ where the coupling with the _L-mode 
is neglected. We conclude that the analysis of the axion- 
photon resonance given in Sect. IIIIBl is accurate. 

We can similarly show that the vacuum resonance dis- 
cussed in Sect. IIII Al is hardly affected by the photon- 
axion coupling. 



D. Numerical Examples of Three-Mode Evolution 

Here we consider the numerical integration of the 
axion-photon system in a neutron star atmosphere, with 
a constant magnetic field B but varying densities. In our 
examples, we assume p = pa exp(— z/i/), where pa is 
the density of the axion-photon resonance, and H is the 
density scale height along the direction of propagation. 
Equation (jSQ]) is equivalent to 



_d_ 
dV 



(72) 



with G given by ((58|) . We integrate (ffS]) outward from 
the high-density region to the low-density region, travers- 
ing the axion-photon resonance at 1^ = 1 and the vac- 
uum resonance at V = 1 + m/q. Fig. [TO] and Fig. [11] 
give two examples of the evolution, with the photon en- 
ergy w = 1 keV (Fig. [TO]) and w = 0.3 keV (Fig. 
respectively, both with the parameters nia = 10^^ eV, 
g = 10-6 GeV"\ B = 10" G, H = 5 cm, and 
9 — 45°, and the same initial conditions a,t V = 1.4: 
E\\ = 0.9, Ei_ = V0J9 and a = 0. The lower panels 
of Figs. [TOKTTl depict \Ai\^ {i = 1,2,3), where the mode 
amplitude Ai is defined by 



(73) 



10 




0.8 



0.6 



0.4 



0.2 



B 0.8 
< 



0.6 
0.4 
0.2 





V=v /q 



0.5 1 

V = v /q 



1.5 



FIG. 10: The evolution of the three-state system in an in- 
homogeneous neutron star atmosphere. The upper panel 
shows the square amplitudes \a\'^, and the lower 

panel shows the square amplitude of non-crossing eigenmodes, 
1^1 1 ^, 1^2 |^3p- The evolution starts in the high-density 
(large Ve — uj'pe/uP) region, where amplitudes are set to 
= 0.9, E±_ — \/0.19 and a — Q. The parameters are 
rua = 10-^ eV, g = 10"'' GeW-\ B = lO" G, ii" = 5 cm, 
61 = 45°, and = 1 keV. 



and $i are the normalized mode eigenvectors {'^^'^i — 1) 
discussed in Sect. IIII Cl and Ai = <i>^$. 

The numerical results fully agree with the analytical 
consideration of Sect. IIII C) and show that the two res- 
onances can be treated independently of each other. In 
particular, the mode evolution across the resonance is de- 
scribed by the non-adiabatic jump probability to within 
about 30%, Eq. ([55]) (for vacuum resonance) or (for 
axion-photon resonance) for values of 7rcs that vary by 
nearly seven orders of magnitude from 10^^ to 10. 



IV. AXION-PHOTON PROPAGATION IN THE 
SUN AND WHITE DWARFS 

The central core of the Sun may a source of keV axions. 
These axions are the target of the CAST experiment Q . 
Could the axions oscillate/transform into photons as they 
travel across the solar envelope and atmosphere? 

The general axion-photon propagation equation (j39p 
can be applied here, and we have Me <C 1 and A|| is 
negligible. As in Sect. IIII Bl it is adequate to consider 



Eq. 



with A, 



pe 



/(2w). The axion- 



FIG. 11: Same as Fig. 10. except uj = 0.3 keV. 

photon mixing angle is given by [see Eq. (|46p ] 

Am 



tan20„i — 



(A,- Ap)/2- 



(74) 



Since Am/ I A 



5.92 xlO~ 



r^^gWiBi sin 9 (where Bi 
is the field strength in units of Gauss), we see that the 
mixing angle is quite small except near the resonance, 
where Ap = Aa- The resonance occurs when = Wpo, 
at the density p ~ 1.21 x \Q~'^Y~^m\ g cm~^. At the 
resonance, the adiabaticity parameter is 



4A2,iJ , Am a 



(75) 



where Rq is the solar radius and H is the density scale 
height. Clearly, 7res ^ 1 and no resonant conversion of 
axion to photon is expected around the resonance. 

In the case of magnetic white dwarfs, with B in the 
range of 10^-10^ G, we still consider <C 1. The axion- 
photon resonance is at Ve = gsin^6' -I- {nia/toY, or the 
density 

p,cs = 1.21 X IQ-^Y-^ml (1 + A||/|A„|) g cm'^, (76) 

where A||/|Aa| = 1.855 x 10^ tj^B^ml^ sin^ 6*, and Bg = 
5/(10^ G). The adiabaticity parameter at the resonance 
is 



7r. 



|A„|(1 + |A||/|A,|) 



2.49 X 10^ 



LJi (m^ ^ggSg sin / H 



|A„ 



10-2i?, 



■0 - 



(77) 
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where we have scaled the density scale height H to the 
typical radius of the star, lO~^i?0. Thus adiabatic reso- 
nant axion-photon conversion is a possibility. For a given 
set of axion parameters, we expect that the flux of pho- 
tons satisfying 7rGs 3> 1 is reduced by a factor of 2, and 
significant linear polarization will be present. Observa- 
tions of the spectra and polarizations of magnetic white 
dwarfs would provide constraint on the axion parame- 
ters. To carry out such constraint quantitatively, it is 
necessary to model the low-density region [see Eq. (1761) ] 
of the white dwarf atmosphere, a task beyond the scope 
of our paper. 

The evolution of the axion-photon system in the vac- 
uum region outside a white dwarf with a varying mag- 
netic field is similar to the neutron star case studied in 
Sect. El 

V. DISCUSSION 

In this paper we have studied the propagation of the 
axion-photon system in the atmospheres and the near 
vicinity of magnetic stars. Because of their strong mag- 
netic fields, we have focused on neutron stars, but our 
main results/methods can be similarly applied to other 
astrophysical bodies. Our study goes beyond previous 
work |11| in that we quantify the various parameter 
regimes for which axion-photon conversion is important, 
we calculate the conversion probabilities for propagation 
in varying magnetic fields and varying plasma densities, 
and we present examples to illustrate how axion-photon 
coupling may affect observed photon spectra and polar- 
izations. 

The axion-photon resonance (where maximum mix- 
ing occurs) is always present in magnetized neutron 
star atmospheres/magnetospheres. This resonance [see 
Eq. (|43p ] is at a higher density than the previously- 
studied vacuum resonance. Complete axion-photon con- 



version is possible only when the adiabatic condition 
[7rcs ^ 1; see Eq. ([5^ ] is satisfied. A high axion-photon 
coupling strength, gradual density gradient and lower 
photon energy tend to make such conversion a possibil- 
ity. Note that in this paper we have treated the plasma 
dielectric tensor in the cold-plasma approximation. For 
low-energy photons (e.g. optical), the resonance occurs 
in the neutron star's magnetosphere and the cold plasma 
treatment is no longer valid. It would be straightforward 
to to generalize our results to more general plasma di- 
electric tensors. 

Even without the axion-photon resonance, partial 
axion-photon conversion may take place during the prop- 
agation in the vacuum region with spatially varying mag- 
netic fields (see Sect. II). 

Applying our result to the axions produced at the cen- 
ter of the Sun, we find that there is no possibility for 
appreciable axion-photon conversion during propagation 
(see Sect. IV). Thus the nondetection of axions in the 
CAST experiment cannot be explained by such oscilla- 
tion effect. Our analysis also shows that with the ax- 
ion parameters allowed by the PVLAS experiment, sig- 
nificant photon-axion resonant conversion is possible in 
highly magnetized white dwarfs. This may produce in- 
teresting spectral and polarization signatures in the ob- 
served radiation from the white dwarf. 
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